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Electromagnetic energy loss of tau leptons is an important ingredient for eventual tau neutrino 
detection from high energy astrophysical sources. Proposals have been made to use mountains 
as neutrino converters, in which the emerging tau decays in an air shower. We use a stochastic 
evaluation of both tau neutrino conversion to taus and of tau electromagnetic energy loss. We 
examine the effects of the propagation for mono-energetic incident tau neutrinos as well as for several 
neutrino power-law spectra. Our main result is a parameterization of the tau electromagnetic energy 
loss parameter /3. We compare the results from the analytic expression for the tau flux using this j3 
with other parameterizations of /?. 



I. INTRODUCTION 

The observed neutrino fluxes in solar neutrino exper- 
iments 0] and in atmospheric neutrino experiments Q 
lead to a picture of neutrino oscillations with substantial 
— > v T mixing 3]. For high energies over astronomical 
distances, astrophysical sources of v e and can become 
sources of v T as well. One method of detecting v T from 
astrophysical sources is to look for the charged particles 
they produced in or near neutrino telescopes. 

Because of the short lifetime of the tau, v T — > r — * v T 
regeneration can be an important effect as the v T passes 
through a significant column depth through the Earth 
H IE IE 0- 0- This is typically relevant for neutrino 
energies below E„ ~ 10 8 GeV. At higher energies, one 
can i mag ine usin g th e Earth as a tau neutrino converter 
II SEfl IH 13 111 El 13 El III 13 III Because of 
the relatively short interaction length of the neutrino at 
high energies, column depths on the order of the distance 
through a mountain or from small skimming angles rel- 
ative to the horizon are most important. Following the 
tau neutrino conversion, a tau is detected by its decay 
e.g., in ice with IceCube, RICE or the proposed ANITA 
detectors, or after it emerges into the air. The possibil- 
ity of detectable signals depends on the incident fluxes, 
of which there are a variety of models [2(|, and on the 
detectors and their sensitivity to various energy regimes. 

Theoretical calculations often rely on approximating 
the neutrino interaction and tau energy loss in terms of 
analytic functions of energy and column depth [Til Il2| . 
These analytic tools are important in exploring the pos- 
sibility for tau neutrino signals in a variety of detectors, 
even though ultimately a full simulation of tau produc- 
tion and energy loss is required. In the energy regime 
where E v ^> 10 s GeV Huang, Tseng and Lin have made 
an evaluation of the smearing of energy due to tau prop- 
agation by doing an approximate stochastic evaluation of 
electromagnetic energy loss of the tau El- Looking at 
fixed incident neutrino energies, they find that the rel- 



ative energy fluctuation of the tau can become large at 
high energies, especially for larger column depths. This 
has the potential to make interpretation of fluxes difficult 
at high energies. 

In this paper, we evaluate tau neutrino and tau prop- 
agation using a stochastic evaluation of tau electromag- 
netic energy loss based on the work of Ref . [3] . We find 
the approximate evaluation of Ref. El reliable. In addi- 
tion to considering mono-energetic incident neutrinos, we 
also consider incident neutrino fluxes with a power law 
behavior, E~ n . We compare the flux of taus following 
propagation through 10 km rock and 100 km rock using 
the stochastic evaluation with estimates of the emerging 
flux of taus using analytic approximations. We refine 
some recent parameterizations of the tau energy loss as 
a function of tau energyjllj, [12( to better approximate 
the stochastic result for the flux of emerging taus using 
analytic approximations. The best parameterization of 
the electromagnetic energy loss parameter /3 is with a 
logarithmic dependence on tau energy. 

In the next section, we describe the models for the high 
energy neutrino nucleon cross section and the tau energy 
loss used in our Monte Carlo simulation. Results of the 
simulation and analytic approximations for the tau en- 
ergy as a function of distance traveled, for monoenergetic 
neutrinos, appear in Section III. Section IV includes re- 
sults and analytic approximations for power law incident 
neutrino spectra. A summary of our results and a com- 
parison with a simplified analytic approach are in the 
final section. 



II. NEUTRINO INTERACTIONS AND TAU 
ENERGY LOSS 

Our evaluation of neutrino interactions and tau en- 
ergy loss is performed using a one-dimensional Monte 
Carlo computer simulation. The neutrino charged cur- 
rent and neutral current interactions are calculated based 
on neutrino- isoscalar nucleon cross sections [3 [3] using 
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FIG. 1: Neutrino interaction length L- lnii in rock as a function 
of neutrino energy. 



is evaluated similarly. We indicate the energy of the 
tau at the point of its production by E\. The neu- 
trino interaction probability and resulting tau energy is 
evaluated using a one-dimensional Monte Carlo program. 
In our analytic approximations described below, we use 
(y) ~ 0.2, which is approximately correct at ultra-high 
energies [%| . An indication of this is shown by the solid 
line in Fig. [3 which shows the ratio of the average tau 
energy to the initial neutrino energy, as a function of 
neutrino energy, in the absence of electromagnetic energy 
loss corrections. 

Our evaluation of the tau electroma gne tic energy loss 
follows the procedure described in Ref . 1211 . Based on the 
scheme outlined by Lipari and Stanev |28| and others , 
we have incorporated the electromagnetic mechanisms of 
ionization, bremsstrahlung (brem), e + e~ pair production 
(pair) and photonuclear (nuc) scattering. The tau life- 
time is taken into account. 

The average energy lost per unit distance is commonly 
described by the formula 



the CTEQ6 parton distribution functions (PDFs) 24]. 

At ultra-high energies, there is a competing effect in 
the neutrino-nucleon cross section as a function of Q 2 = 
~q 2 > 0, the quantity describing the four-momentum 
squared of the vector boson. PDFs increase as a func- 
tion of InQ 2 (2{|, however, the vector boson propagator 
effectively cuts off the growth in Q 2 at the vector boson 
mass. For charged current scattering, this means that the 
value of parton x = Q 2 /2p-q for nucleon four-momentum 
p is approximately 



ME,, 



(1) 



in terms of the incident neutrino energy E v , nucleon mass 
M and W boson mass My/. For E v — 10 12 GeV, one 
finds that x ~ 10~ 8 . 



The CTEQ6 PDFs are parameterized for 10" 



< x < 



1, so one must extrapolate for smaller values of x. The 
PDFs at small values of parton x (x < 10 -6 ) are extrap- 



olated here using a power law xq(x, Q ) 



to DGLAP evolved PDFs at x„ 



10 



-6. 



xq(x 



matched 



(2) 



This approach is based on the result [26| that for A ~ 
0.3 — 0.4, the gluon PDF has the approximate form of 
xg(x,Q) ~ x~ x , and g — > qq splitting is responsible for 
the sea quark distributions that dominate the cross sec- 
tion at ultra-high energies. We show in Fig. ^the neu- 
trino interaction length in rock. Other approaches to 
the small x extrapolations for ultra-high energy neutrino 
cross sections yield similar interaction lengths to 
within a factor of approximately 2 ±:L at 10 12 GeV. 

The differential distribution for neutrino nucleon scat- 
tering in terms of the lepton inelasticity y, 



V 



E v — E\ 
E v 



(3) 



(^) = ~(a + pE T )p 



(4) 



In this expression, z is the distance the tau travels. The 
quantity a ~ 2 x 10~ 3 GeV cm 2 /g accounts for ionization 
energy loss [3(3]. Here, (3 is 



13 = f3 Diem + /3 pair + /3 nu ' 
^N A f da' 



(5) 



for Avogadro's number Na, atomic mass number A, tau 
inelasticity y and i = brem, pair and nuc. 

In our treatment here, as in Ref. we separate (3 
into continuous and stochastic terms: 



P — Pcont + Pa 



N A 



A 



10" 



da 1 N A 



= — r / d yy-r- + — r 



dy 



A 



dyy-r- ■ (6) 



io~ 3 



dy 



We treat the tau propagation stochastically, with a se- 
ries of electromagnetic interactions occurring according 
to probabilities based on da 1 /dy for y > 10 -3 and the de- 
cay probability. For each distance between interactions in 
the Monte Carlo program, the tau energy loss calculated 
from Eq. (3) is supplemented by Eq. Q with the substi- 
tution of j3 — > Pcont ■ By doing it this way, we account for 
the full electromagnetic energy loss, but we limit com- 
puter time spent where the cross section is large, but the 
energy loss is small (y < 10~ 3 ). The detailed formu- 
las for the bremsstrahlung [3l|], pair production and 
photonuclear [33| processes are discussed and collected in 
Ref. |2l|. See also Ref. [34[ for the muon case. 

For the muon, [3 is dominated by pair production and 
bremsstrahlung for E^ = 10 3 — 10 7 GeV. At higher ener- 
gies, (3 grows slowly with energy due to the photonuclear 
process which contributes a larger fraction for E^ > 10' 
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GeV. For the tau, the photonuclear process is compa- 
rable to pair production, while the bremsstrahhmg con- 
tribution is suppressed, leading to an increase in (3 with 
energy for E T > 10 3 GeV. 

In the analytic approximations discussed below, we will 
consider several forms of as a function of E T , and we 
will treat 



[ dz ' 



dz 



(7) 



To guide our discussion, it is useful to compare the rele- 
vant constants and distance scales. First, at 10 8 GeV, 



(3E T ~ 60 



cm 2 GeV 



10 8 GeV 



(8) 



is much larger than a, so we can essentially ignore a. 
The tau range is then characterized roughly by 



1 



6 km 



(9) 



in standard rock (p = 2.65 g/cm 3 ). In fact, the range is 
somewhat larger than this estimate. The decay length of 
the tau is 



7cr ~ 5 km 



10 8 GeV 



(10) 
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FIG. 2: Average tau energy scaled by the energy of incident 
neutrinos that produced the taus, without energy loss (solid), 
and with energy loss for taus produced in 10 km and 100 km 
rock (dashed). 
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It is in the energy regime near 10 8 GeV that the transi- 
tion from the lifetime dominated to the energy loss dom- 
inated range for the tau occurs. 



III. MONO-ENERGETIC NEUTRINOS 
A. Results from Monte Carlo simulation 

We begin our discussion by considering single incident 
neutrino energies to see the effect of energy loss com- 
bined with decay probabilities for two different trajecto- 
ries: through 10 km rock (26.5 km.w.e.) and through 100 
km rock. Fig. [21 shows the average tau energy divided 
by the incident neutrino energy as a function of incident 
neutrino energy. The average distance that the tau trav- 
els after production which survives to emerge from the 
rock (10, 50 and 100 km) is shown in Fig. [3J 

For 10 km of rock, above 10 9 GeV, the average dis- 
tance the tau travels is 5 km. The neutrino interaction 
probability is equal anywhere in the 10 km path, and the 
energies are large enough that the decay length is long 
compared to the path. For 50 km and 100 km, Fig. |3| 
shows that the last ~ 15 km are the most important for 
the emerging taus. 

The standard deviation, normalized by the average tau 
energy, 



V{E 2 T ) {E T f 
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FIG. 3: Average distance the tau travels after production 
which survives to emerge from 10 km rock (solid), 50 km 
rock (dashed) and 100 km rock (dot-dashed), as a function of 
incident neutrino energy. 



is shown is Fig. Higher energies and larger distances 
lead to larger fluctuations. Huang et al.'s results in Ref. 
[T^| correspond to our Figs. [21 and 0] and are in good 
agreement. We comment that by normalizing the stan- 
dard deviation to the incident neutrino energy rather 
than average energy of the emerging tau, one finds a more 
stable result for the ratio. The standard deviation scales 
approximately as (0.20 — 0.25) x E v for 10 km rock over 
the same energy range as the figure, and with a factor of 
(0.13 - 0.25) x E v for 100 km rock. 

A final result for mono-energetic neutrinos is shown in 
Fig. [SI Each histogram shows the average energy of the 
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FIG. 4: Standard deviation a scaled by the average tau energy 
as a function of the incident neutrino energy without energy 
loss (solid) and for 10 km and 100 km of rock (dashed). 
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The solution depends on the energy dependence of 0. 
Several parameterizations have been used, including a 
linear dependence on energy [l2T | and a logarithmic en- 
ergy dependence 8]. We consider three cases: 



(I) 



(II) E T 



E T = E l r e-PP z ' 
= constant 



-@TPZ f 



/3 x +7 T B*(l-e-^f 2 
= T + lT E T 



(13) 



(III) E T = exp 



/3i 



D -/3ipz' 



+ \n(Ei/E ) 



E 



= 00+13! HE/E ) . 

The quantity z' is the distance traveled by the tau af- 
ter its production. None of the choices for is entirely 
satisfactory. Case (III) is the best approximation to the 
energy curves of Fig. HJ1 For case (I), 



= 0.85 x 1CT 6 cm 2 /g 



(14) 



works moderately well for E v = 10 10 GeV, but poorly for 
higher and lower energies. We reproduce the histograms 
of the stochastic results from Fig. showing the last 20 
km in Fig. Altogether with parameterizations of 0. In all 
cases of the parameterized energy, we take E\ — 0.8£?„. 
The dot-dashed lines show the parameterization of case 
(I)- 

Case (II) is used by Aramo et al. in Ref. with 



T = 0.71 x 10~ 6 cm 2 /g 

7r = 0.35 x 10~ 18 cm 2 /g GeV 



(15) 



FIG. 5: Average energy of emerging tau as a function of its 
production point for a total distance of 100 km of rock, for 
fixed initial neutrino energies between 10 s — 10 12 GeV. 



emerging tau as a function of the production point z for a 
total depth of 100 km and at fixed initial neutrino energy. 
The energies of the incident neutrinos are E v — 10™ GeV 
for n = 8, 9, 12 (lower to upper curves). 



B. Analytic approximations 

Standard approaches to parameterizing the final en- 
ergy in terms of the initial energy begin with 

^ ^ -0pE T (12) 



The energy parameterization in this case is shown with 
the dashed lines in Fig. 

The parameterization using a logarithmic dependence 
on energy for does the best at reproducing the energy 
behavior of the emerging taus. This is shown with solid 
lines in Fig. with 

00 = 1.2 x 10~ 6 cm 2 /g (16) 

01 = 0.16 x 10~ 6 cm 2 /g 
E = 10 10 GeV . 

Not surprisingly, this energy dependent will give the 
best representation of the Monte Carlo simulation results 
for the emerging tau fluxes from power law incident neu- 
trino fluxes. We note that these values of are larger 
than what one would extract from the direct calculation 
of [2l| as in Eq. (6). This is because Eq. (7) is only 
an approximate equality. 
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FIG. 6: Average energy of emerging tau as a function of its 
production point for a total distance of 100 km of rock, for 
fixed initial neutrino energies between 10 s - 10 12 GeV from 
Monte Carlo simulation (histograms), and parameterization 
of (3 as in case (I) (dot-dashed), case (II) (dash) and case (III) 
(solid). 



FIG. 7: Flux ratio of tau flux to incident neutrino flux as 
a function of energy for 10 km and 100 km of rock, for an 
incident neutrino energy scaling like E^ 1 . The dashed lines 
do not have energy loss or decay included. The histograms 
include the stochastic treatment of the neutrino interaction 
and tau energy loss and decay. 



IV. POWER LAW NEUTRINO SPECTRA 
A. Results from Monte Carlo simulation 

In this section, we consider cases of incident neutrino 
fluxes that depend on energy. Rather than use specific 
flux models, we consider power law fluxes F v ~ E~ n , for 
n = 1,2,3, to see qualitatively and quantitatively how 
flux propagation is affected by tau energy loss. 

In Fig. [7| we show the results of the stochastic prop- 
agation of the incident tau neutrino, which converts to 
a tau and continues through the depth D of rock, for 
D = 10 km and 100 km and an incident neutrino flux 
scaling like E^ 1 . All of the fluxes shown in this section 
have a sharp cutoff at E v = 10 12 GeV. The two upper 
dashed curves show the outgoing tau flux to the incom- 
ing neutrino flux neglecting tau energy loss and decay. 
Because of the larger column depth for 100 km of rock, 
the ratio of the tau to neutrino flux is ten times larger 
than for 10 km. The flux ratio increases with energy to 
scale with the neutrino cross section. 

Our stochastic propagation including energy loss and 
decay is shown by the histograms in Fig. Energy 
loss is much stronger for the 100 km depth, so the ratio 
of outgoing tau flux to incident neutrino flux is smaller 
for 100 km than 10 km for E > 10 9 GeV, despite the 
increased number of targets for the neutrinos. A quali- 
tative discussion of the results for Fig. [7] is in Sec. IV. B. 

Figs. 03 and show the ratio of the outgoing tau flux 
to the incident tau neutrino flux for incident neutrino 
fluxes scaling as E~ 2 and E~ 3 , respectively. The flux 
ratios including decay and energy loss at E — 10 8 GeV 
for these steeper fluxes can be understood based only on 



the lifetime, since energy loss is not overwhelming at that 
energy and there is little pile-up of taus from neutrino 
interactions at higher energies, followed by tau energy 
loss. The dashed lines have a factor of D to account 
for the column depth of nucleon targets in the neutrino 
interaction probability. When decays are included, the 
relevant distance at E = 10 8 GeV is d = E t ct/tti t , so 
the neutrino interaction probability is proportional to d 
rather than D. Neutrino attenuation in 10 km and 100 
km of rock at E v ~ 10 8 GeV is small, so the flux ratios in 
Figs. |Sland El are the same for the two depths. At higher- 
energies, attenuation and energy loss play a bigger role. 

Our aim here is to compare the stochastic propagation 
to analytic approaches. In the next section we outline the 
analytic evaluation of the emerging tau flux and make 
quantitative comparisons with the histograms in Fig. \7\ 
E and M 



B. Analytic approximations 

The analytic approximation to the emerging tau flux 
as a function of tau energy depends on a number of com- 
ponents. We review the discussion of Ref. [TlJ here, ex- 
panded to include two energy dependences of 0, a linear 
dependence [l^ and a logarithmic one Q . 

The incident flux F V (E V , 0) is attenuated along its tra- 
jectory a distance z through the rock with an approxi- 
mate factor 



F u (Eu,z) ^exp -za tot (E v ) P N A F u (E v ,0) 



(17) 



In our evaluation here, we are assuming a constant den- 
sity (standard rock). Neutrino regeneration0,0,0,0,l3- 
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FIG. 8: Flux ratio of tau flux to incident neutrino flux as a 
function of energy for 10 km and 100 km of rock, for an inci- 
dent neutrino energy scaling like E~ 2 . The dashed histograms 
do not have energy loss or decay included. The histograms 
include the stochastic treatment of the neutrino interaction 
and tau energy loss and decay. 
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FIG. 9: Flux ratio of tau flux to incident neutrino flux as 
a function of energy for 10 km and 100 km of rock, for an 
incident neutrino energy scaling like E^ 3 . The dashed lines 
do not have energy loss or decay included. The histograms 
include the stochastic treatment of the neutrino interaction 
and tau energy loss and decay. 



from neutral current processes or tau decay is not impor- 
tant over these short distances. 

The probability for neutrino conversion to taus in the 
distance interval [z, z + dz] with energy in the interval 
[E\,E\ + dE* T ] is 



Pcc{vt -» r, E v ,Ei) = dzdE l T N A p 



do cc (E Vl E\) 
dEl 



(18) 



As discussed above, we take E\ = 0.8^ 
dace 



dE\ 

(y) 



a CC {E v )8{Ei - (1 - {y))E v ) (19) 
0.2 . (20) 



Writing the differential cross section this way ensures 
that the integral over tau energy yields the charged- 
current (CC) cross section. 

In the infinite lifetime, no energy loss (r — > oo, [3 — > 0) 
limit, the emerging tau flux given incident tau neutrinos 
is 



F T (-Ey)i"— > 



oo,/3^0 







x exp 



dz 



ii 



zG tot {E v )pN A \F v {E v ^) 

dF?(E* T1 z) 
dz 



(21) 



for rock depth D. 

In the case where the energy loss is treated analytically, 
the lifetime and energy loss are accounted by a factor 
P S urv{E l T ,E Tl D — z) as the tau travels over the remaining 
distance from production point z to emerge from the rock 
after distance D with outgoing energy E T . In fact, E\ : 
E T and D — z are not all independent, as shown in Eq. 
(|13fl . The survival probability as a function of energy 
comes from the solution to the equation 



dP„ 



dz 



ctE t /t 



(22) 



together with the approximate dE T /dz formula of Eq. 
O , leading to 



dP, 



dE T 



m r p 

crPpE* surv 



(23) 



The solution depends on the energy dependence of (3. For 
constant (3, the solution is 



P S urv{E T ,E l T ) = exp 



m T 



1 



cr(3p \ E T 



The emerging tau flux is 

F T {E T ) = [ 
Jo 



case (I) . 
(24) 



dz 



Psurv (-^T ) ) 



dEl 



x 5(E? - El (E T , D - z) J ~^dK 

D dF°(E 1 ' z) 
dz 







x 6(E T -E T (Ei!,D-z))dE% 



(25) 



The delta function explicitly enforces Eq. I|13|l rewritten 
with E\ in terms of E T and D — z. A factor of dE"/dE T 
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FIG. 10: Flux ratio of tau flux to incident neutrino flux as a 
function of energy for 10 km of rock, for an incident neutrino 
energy scaling like E^ 1 . The histograms include the stochas- 
tic treatment of the neutrino interaction and tau energy loss 
and decay. The analytic approximations for cases (I, dot- 
dashed), (II, dot-dot-dot-dashed) and (III, dashed) are also 
shown. 



FIG. 11: Flux ratio of tau flux to incident neutrino flux as a 
function of energy for 100 km of rock, for an incident neutrino 
energy scaling like E~ 2 . The histograms include the stochas- 
tic treatment of the neutrino interaction and tau energy loss 
and decay. The analytic approximations for cases (I, dot- 
dashed), (II, dot-dot-dot-dashed) and (III, dashed) are also 
shown. 



accounts for the fact that F T (E T ) represents the number 
of r's per tau energy interval, say for AE T = E% — E\. 
Electromagnetic energy loss means that one must sample 
a much larger interval of initial tau energies, AE\ = 
E2 exp(f3pD) — Ei. Combining this factor with the delta 
function leads to the factor of 6(E T — E" exp(—/3p(D — 
z))) that appears in the literature|ll| for tau neutrino 
production of taus for a constant beta: 



F T {E T ) = [ dz 
Jo 



dF°(E l T ',z) 



dz 



~ Psurv i^E T i E T ) 

6(e t - E% exp(-(3p(D - z)))dE* . (26) 



In the f3 — > limit, the emerging tau flux is the ex- 
pected result, 



F T (E T ) 



(i 



D dF\E^z) , 

dz exp 

dz v 



m T {D 



ctE t 



(27) 



Continuing with case (I) for non-zero (3, explicit sub- 
stitutions give 

F T {E T ) = J dz J dE v exp(-zatot(E„)pN A )F„(E 1/ ,0) 
x N A pacc{E v )8{E T - 0.8£„ exp(-/3p(D - z))) 



x exp 



crf3p \ E T 



(28) 



A useful limit to understand our results in Fig. \7\ is 
the case of constant j3 in the long lifetime limit. For 
simplicity, we set (y) = so that E l T — E u 11]. If we also 



assume D <C Lint where L- lnt is the neutrino interaction 
length, we can simplify Eq. to 



F T (E T ) = / / dE v F v (E v , 0)5(E T - Eve- 

Jo L iat J 



0p(D-z) 



E T ppLgcj E . 



dE v F v {E v ,Q), 



(29) 



for Li n t approximately constant. The energy integral 
runs from E m i n — E T to E max = min(E T exp(/3pZ?), 10 12 
GeV). At low energies and short distances, (e.g., E T ~ 
10 s GeV and D = 10 km), the maximum energy is 
E T exp(ppD). For F u ~ E~\ F T ~ {D/Lgf )-E-\ while 
for F v - E~ n for n > 1, F T ~ l/((n - l)PpLg° ) ■ E~ n . 
This means that even with strong energy loss effects, for 
an E^ 1 flux in the appropriate energy regime, the total 
distance D determines the tau flux rather than the char- 
acteristic scale of the energy loss, l//3p. In the actual 
case of the r, the tau decay length jct is comparable to 
the energy loss scale l/(3p at E = 10 8 GeV. For the E^ 1 
flux, D is still the dominant scale, although the details 
including both energy loss and decay as well as energy 
dependent cross sections have an effect, as seen in Fig. 
H Figs. |H1 and El show that this is not the case for E~ 2 
and E^. 

Eq. (|25|l is the master equation for all three cases 
discussed here for the three parameterizations of (3. The 
argument of the delta function comes from Eq. Ijl3|l for 
the specific cases I-III. The survival probabilities come 
from the solution to Eq. I|23|) for the appropriate energy 



8 



dependent p. This gives (see Ref. [13 for case (II)) 
' E\(p r + lT E T )\^ 



Psurv \Et 7 ) — 



E T {f3 T + ~f T Ei) 



x exp 



1 



1 



Psurv {E T ,E l T ) = exp 



ctI3 tP \ E T El 



case (II) 



E T \E 



E. 



x exp 



1 ( , K 



nij 



CTf3 p \ E T El 



case (III) . 
(30) 



In case (III) , an expansion in Pi / Pq has been made in Eq. 
(|23[1 . Numerically, we have kept terms through Pf/P^, 
however in Eq. 130J1 . we show the result where only the 
first term has been kept. 

None of the parameterizations can completely describe 
the Monte Carlo results. Based on Fig. 6, it is not sur- 
prising that the form of (3 which depends on log E T /E 
works the best to describe the outgoing tau flux. In Figs. 
1101121 we show a comparison of the various parameteriza- 
tions of the outgoing tau flux, normalized to the incident 
neutrino flux. 

Fig. shows the Monte Carlo results for 10 km of 
rock and an incident flux with E~ x energy behavior with 
the histogram and the case (III) parameterization with 
the dashed line. The case (I) parameterization, with con- 
stant (3 appears as a dot-dashed line, while the case (II) 
power law is shown with a dot-dot-dot-dashed line. One 
sees the turnover in the flux ratio as the change is made 
from E max — E T exp(f3pD) to E rnax = 10 12 GeV for the 
E^ 1 flux. The parameterization of f3 as a function of 
\og(E T / Eq) (case (III)) gives results for the flux ratio to 
within about ±50% of the Monte Carlo result. The other 
two parameterizations do not do as well, underestimating 
by a factor of ~ 2 for E = 10 9 GeV and overestimating 
by - 3 for E = 10 11 GeV. 

For 100 km of rock, with the same incident flux, 
the case (III) parameterization compares well (under- 
estimates by ~ 15%) the Monte Carlo result between 
10 9 — 10 11 GeV. By comparison, the other two parame- 
terizations for 100 km rock overestimate the tau flux by 
a factor of 1.3 — 1.6 and 1.8 — 2 at energies 10 10 and 10 11 
GeV respectively. 

For the E~ 2 incident neutrino flux, the logi? param- 
eterization of (3 does well for 100 km as seen in Fig. 
ITTl Again, the dot-dashed and dot-dot-dot-dashed lines 
show cases (I) and (II) respectively. For incident neu- 
trino fluxes scaling with E~ 3 , the 10 km result is well 
represented by the log-E parameterization, as shown in 
Fig. The overestimate of the flux for cases I and 

II occurs because energy loss is underestimated at high 
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FIG. 12: Flux ratio of tau flux to incident neutrino flux as a 
function of energy for 10 km of rock, for an incident neutrino 
energy scaling like E~ 3 . The histograms include the stochas- 
tic treatment of the neutrino interaction and tau energy loss 
and decay. The analytic approximations for cases (I, dot- 
dashed), (II, dot-dot-dot-dashed) and (III, dashed) are also 
shown. 



energies with these parameterizations of p. 



V. DISCUSSION 

We have performed a one-dimensional Monte Carlo 
simulation of v T conversion to r and r propagation in- 
cluding electromagnetic energy loss over distances of 10- 
100 km in rock. Analytic approximations are useful tools 
in exploring the possibility of detecting a variety of pre- 
dictions for incident neutrino fluxes 11]. Using mono- 
energetic sources of tau neutrinos, we modeled the tau 
energy loss with an energy loss parameter (3 depending on 
\og(E T / Eq). In an analytic approximation for the outgo- 
ing tau flux, this model for energy loss does better at rep- 
resenting the Monte Carlo results than two other choices 
for (3: a constant and a (3 increasing linearly with tau en- 
ergy. Our improved parameterization of (3 will help refine 
theoretical predictions of tau neutrino induced events in 
neutrino telescopes. 

It is not difficult to numerically evaluate the integrals 
of Eq. H25|) including the more complicated logarith- 
mic dependence on energy of (3 (case (III)) for P surv 
and E Tl nevertheless, it is interesting to see the impact 
of the modifications from that energy dependence. Al- 
though formally inconsistent, one could simply use the 
flux calculation for constant (3, Eq. H2b|) . and substitute 
p = p + p 1 \n(E T /E ). This is shown in Fig. U3J for 
D = 100 km and E~ , with the dot-dashed line. As 
in Figs. 1101121 the dashed line comes from the full case 
III (logarithmic dependence) result, and the histogram 
represents the Monte Carlo result. The dot-dashed line 
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FIG. 13: Flux ratio of tau flux to incident neutrino flux as a 
function of energy for 100 km of rock, for an incident neutrino 
energy scaling like E^ 1 . The histogram includes the stochas- 
tic treatment of the neutrino interaction and tau energy loss 
and decay. The analytic approximations using Eq. (28) with 
P = Po + /3i ln{E T /E ) (dot-dashed), using E T {E\,D - z) for 
case III but a simplified P 3U rv (triple dot-dashed) and the full 
case III result (dashed) also shown. 



is nearly identical to the complete case III evaluation. 
The triple dot-dashed line uses the full expression for 



the energy dependence on distance from Eq. i|13|) . but 
it uses the simplified form of the survival probability: 
P surv = exp[-m T /(cT/3p) ■ (E- 1 - Ey 1 )]. This does 
as well in representing the Monte Carlo result. Of the 
fluxes and rock depths that we have discussed here, Fig. 
1131 is representative of the three curves incorporating to 
varying degrees the logarithmic energy dependence of (3. 

In summary, our Monte Carlo results for tau neutrino 
interaction and tau propagation through rock suggest 
that the tau electromagnetic energy loss parameter be 
parameterized by /3 = fio+fii \og(E/E ) (case III), where 
the parameters /3q, (3\ and En appear in Eq. (16). Using 
the analytic formula of Ref. [ll| , modified to account for 
the logarithmic energy dependence, a reasonable agree- 
ment with the Monte Carlo results is obtained. Using 
the formulae of Ref. for constant /?, with the small 
change of taking (y) — 0.2 but substituting @(E T ), also 
is in reasonable agreement with the Monte Carlo result. 
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